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Abstract: We address the problem of achieving a random laser with 
a cloud of cold atoms, in which gain and scattering are provided by the 
same atoms. In this system, the elastic scattering cross-section is related 
to the complex atomic polarizability. As a consequence, the random laser 
threshold is expressed as a function of this polarizability, which can be 
fully determined by spectroscopic measurements. We apply this idea to 
experimentally evaluate the threshold of a random laser based on Raman 
gain between non-degenerate Zeeman states and find a critical optical 
thickness on the order of 200, which is within reach of state-of-the-art 
cold-atom experiments. 
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1. Introduction 

Random lasing occurs when the optical feedback due to multiple scattering (or "radiation trap- 
ping") in a gain medium is strong enough so that gain in the sample volume overcomes losses 
through the surface. Since its theoretical prediction by Letokhov HI, great efforts have been 



made to experimentally demonstrate this effect in different kinds of systems J2]|3]|4j|5]|6l, as 
well as to understand the fundamentals of random lasing Q [H |9l [TO) . The broad interest of 
this topic is driven by potential applications (see fTT) and references therein) and by its con- 
nections to the fascinating subject of Anderson localization lfl2l . State-of-the-art random lasers 
ifTTl [131 [T4l are usually based on condensed matter systems, and feedback is provided by a 
disordered scattering medium, while gain is provided by an active material lying in the host 
medium or inside the scatterers. In general, scattering and gain are related to different physical 
entities. 

Another system that can be considered for achieving random lasing is a cold atomic vapor, 
using magneto-optical traps lfT5l . where radiation trapping lfl6l [171 as well as lasing lfl8l [191 
have already been demonstrated. One advantage is the ability to characterize and model the mi- 
croscopic properties of the medium, which can be extremely valuable for a better understanding 
of the physics of random lasers. 

However, in such system, the ability to combine gain and multiple scattering at the same 
time is not obvious, as both should be provided by the same atoms. On the other hand, it has 
been shown recently that the peculiarity of this system leads to a simple condition for random 
lasing in the incoherent regime f 20] . The threshold is indeed defined as a critical on-resonance 
optical thickness bo, which is a function of the complex atomic polarizability a as the single 
parameter l2Tl . This has been used to predict theoretically the threshold of a random laser based 
on Mollow gain, for which the atomic polarizability is analytically known ||2T1|221 . A critical 
bo of the order of 300 has been found. 

In contrast to the ab initio theoretical approach of 12TI . we present here an experimental 
evaluation of the threshold of a random laser. Our method relies on the fact that thanks to 
Kramers-Kronig relations | 23|. the complex atomic polarizability is indeed one single indepen- 
dent parameter, and thus can be fully determined by a spectroscopic measurement. This idea 
is general and could be applied with any gain mechanisms. We demonstrate its usefulness here 
with Raman gain between non-degenerate Zeeman states IT8l [191 l24l [25l . We obtain a critical 
optical thickness on the order of 200, lower than with Mollow gain 12TI . 

2. Measuring the threshold of a random laser with cold atoms 

From Letokhov's diffusive description of light transport in a homogeneous, disordered and 
active medium of size L, we know that the random laser threshold is governed by two charac- 
teristic lengths: the elastic scattering mean free path £ sc (26] [23 and the linear gain length £ g 
(£g < corresponds to absorption or inelastic scattering). In the diffusive regime, defined as 
L ^> 4c, the lasing threshold is reached when the unfolded path length, on the order of L 2 /£ sc , 
becomes larger than the gain length. More precisely, the threshold is given by |[Tl [T3l 



where j3 is a numerical factor that depends on the geometry of the sample (j3 = 1 for a slab, 
j3 = 2 for a sphere, which is the case we consider in the following), and L e ff = rjL is the 
effective length of the sample, taking into account the extrapolation length 11261 . For L > £ sc 
and a sphere geometry, T] = 1+ 2^/ [L/£ sc + 2£] with E, ~ 0.71 GH |29). Note that deeply in 
the diffusive regime (L 3> £ sc ), r) ~ 1 . Another important length scale is the extinction length 
£ ex , as measured by the forward transmission of a beam through the sample, T = e~ L '^ a . The 
extinction length is related to the other lengths by £~ = ^sc' — 

For an atomic vapor, these characteristic lengths can both be computed as a function of the 
atomic polarizability a(co) at frequency (0. The extinction cross-section is indeed given by 
Oex(o) = k x Im[a(o)] and the elastic scattering cross-section by <7 S c(ft>) = k 4 /6k x |a(©)| 2 




(1) 



11301 (k = co/c is the wave vector). Note that the first relation is general to any dielectric medium 
whereas the second one is specific to resonant point-dipole scatterers. The characteristic lengths 
are then £^ sc = p C7 ex sc , where p is the atomic density. The gain cross-section can be defined 
the same way by d~ 1 = p Og. The vapor is supposed homogeneous, as well as the pumping field, 
so that both p and a are position-independent. Even though this is not the precise geometry of 
a cold-atom experiment, it allows us to perform analytical estimations. As we consider only 
quasi-resonant light, we shall use k = ko = C0q/c with (Oq the atomic eigenfrequency. In the 
following, we shall also use a dimensionless atomic polarizability a, defined as a = ax 6n/k^, 
and omit the dependence on CO. We can now rewrite a sc = 0o|6:| 2 and a g = Ob (\oc\ 2 — Im(a)), 
where Ob = dn/k^ is the resonant scattering cross-section, such that the threshold condition, as 
expressed by Eq. (Q]), reduces to 

2n 

pa Q L cf{ = rib Q > .■„,..„ _ , (2) 
y / 3\a\ 2 {\a\ 2 -lm(a)) 

where is the on-resonance optical thickness of the cloud. This condition is valid as soon as 
the medium exhibits gain, i.e. \a\ 2 — Im(a) > 0. 

The threshold condition is thus given by a critical on-resonance optical thickness, which is 
an intrinsic parameter of the cloud, expressed as a function of the complex atomic polarizability 
only, which depends on the pumping parameters. Although the initial condition of Eq. ([TJ in- 
volves two characteristic lengths, we emphasize here that this is really one single independent 
parameter, as real and imaginary parts of the atomic polarizability are related via Kramers- 
Kronig relations 1231 . This point is due to the originality of the system that we are considering, 
in which the same atoms are used to amplify and scatter light. 

This property has two important practical consequences. The first one is that we cannot adjust 
one quantity (for example gain) independently of the other (scattering rate, or vice versa), so 
that the existence of reasonable conditions for random lasing is not obvious. This issue has been 
positively answered recently and it has been shown that random lasing can even occur with a 
low amount of scattering (21) . The second one is that only one quantity has to be measured 
to determine the threshold, as soon as we can measure it for every CO, since Kramers-Kronig 
relations involve integrals over CO. A weak probe transmission spectrum, which we can rewrite, 
with our notations, T(co) = exp [—bo x Im(a(a)))], contains therefore enough information to 
fully characterize a(co) and then to deduce the critical optical thickness. In the following, we 
use this idea with Raman gain. 

Note that without this possibility, measuring independently the two characteristic lengths is 
difficult. Besides the transmission spectrum, one needs another measurement, which can be 
provided by the fluorescence. Nevertheless, the probe fluorescence is small compared to the 
pump one, and inelastic scattering is not easily distinguished from elastic scattering. Despite 
these difficulties, preliminary measurements, in a limited range of parameters, have qualita- 
tively validated the approached based on Kramers-Kronig relations |3~T1 . 

3. Application to Raman gain 

Our experiment uses a cloud of cold 85 Rb atoms confined in a vapor-loaded magneto-optical 
trap (MOT) lfl5l produced by six large independent trapping beams, allowing the trapping 
of a few 10 9 atoms at a density of 10 10 atoms/cm 3 , corresponding to an on-resonance opti- 
cal thickness of about 10. To add gain to our system, we use a pump beam, which is tuned 
near the F = 3 — > F' = 4 cycling transition of the D2 line of 85 Rb (frequency C0q, wavelength 
A = 780 nm, natural linewidth r/27T = 6.1 MHz), with a detuning A = (Op — COq, which can 
be changed via an acousto-optic modulator in a double-pass configuration. The pump beam 
has a linear polarization and a waist larger than the MOT size (a few millimeters) to ensures 
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Fig. 1. (a) Principle of the experiment. We send a weak probe beam on the magneto-optical 
trap (MOT) and the transmission is recorded on a photodetector (PD). The probe frequency 
co is ramped during the acquisition in order to record a spectrum. Another, stronger beam of 
frequency cop is used as a pump, (b) Principle of the Raman mechanism (depicted here for 
a F = 1 ^ F' = 2 transition), (c) Experimental transmission spectra, plotted as a function 
of the pump-probe detuning 8. Without pumping, spectrum (1) shows only the atomic 
absorption. A pump beam of detuning A = — 3.8r and intensity 13 mW/cm 2 , corresponding 
to a Rabi frequency Q. = 2.5T, is added to obtain spectrum (2), which then exhibits a Raman 
resonance in the vicinity of S = 0. The atomic absorption is shifted due to the pump-induced 
light shift and the absorption is reduced due to saturation. 



homogeneous pumping. An additional, orthogonally polarized beam is used as a weak probe 
to measure transmission spectra with a propagation axis making an angle with the pump-beam 
axis of about 17° [Fig.[TJa)]. This small angle, together with the low temperature of our sample 
(~ 100 jJ.K) allows us to neglect any relative Doppler broadening (~ 40 kHz). The probe fre- 
quency co can be swept around the pump frequency with a detuning 8 = co — ©p. Both lasers, 
pump and probe, are obtained by injection-locking of semiconductor lasers from a common 
master laser, which allows to resolve narrow spectral features (this has been checked for earlier 
experiments |[T9l down to 10 kHz). All our experiments are time-pulsed with a cycling time of 
30 ms. The trapping period lasts 29 ms, followed by a dark period of 1 ms, when the MOT trap- 
ping beams and magnetic field are switched off. In order to avoid optical pumping into the dark 
hyperfine F = 2 ground state, a repumping laser is kept on all time. Pump-probe spectroscopy 
is performed during the dark phase, short enough to avoid expansion of the atomic cloud. Data 
acquisitions are the result of an average of 300 cycles. 

Raman gain relies on the pump-induced population inversion among the different light- 
shifted mp Zeeman sublevels of the F = 3 hyperfine level ll24l [25), as depicted in Fig. 01b). 
The optical pumping induced by the n -polarized pump laser leads to a symmetric distribution 
of population with respect to the ntF = sublevel of the ground state, with this sublevel being 
the most populated and also the most shifted, due to a larger Clebsch-Gordan coefficient l32l . 
Atoms are probed with a n -polarized (with perpendicular direction) probe beam, thus inducing 



Am/r = ±1 Raman transitions. Depending on the sign of the pump-probe detuning 8, the pop- 
ulation imbalance induces gain or absorption. Each pair of neighboring sublevels contributes 
with a relative weight depending on the population inversion. In practice however, the levels 
are not separated enough to be well resolved, and only two structures (with opposite signs) 
are visible, one corresponding to amplification for 8 = — 8r and one to absorption for 8 = Sr. 
Note that this situation corresponds to a red detuning for the pump (A < 0) and that the signs 
are inverted for blue-detuning (A > 0). As 8r comes from a differential light-shift (because of 
different Clebsch-Gordan coefficients), it is usually on the order of T/10, whereas A is a few 
r. The width y of the resonances is related to the elastic scattering rate, also much smaller than 
r ll24l . Far from the main atomic absorption resonance, the Raman resonance is thus a narrow 
spectral feature, as in Fig. [TJc), such that we can fit it independently of the main absorption 
line [Fig. [2a)]. Therefore, we scan the frequency of the probe beam around 8 = only, which 
reduces the interaction time with the pump, thus suppressing radiation pressure and subsequent 
unwanted Doppler shift. Note that adding a second counterpropagating pump beam is not a 
suitable solution, as in this situation, other mechanisms may occur (recoil-induce resonances, 
four- wave mixing, atom localization in potential wells) l32l[33l[34l . which would complicate 
the analysis. We use the polarizability 6:r(5,A,£2) to describe the Raman structure, with 



This function is particularly convenient as the Kramers-Kronig transformation of a Lorentzian 
profile is well known. We thus avoid any numerical integration. 

Our experimental procedure is the following. We scan the probe frequency from 8 = — F to 
8 = r during 100 /is and record one Raman transmission spectrum. During the same cycle, we 
perform two larger scans without pumping, one before the pump-probe spectroscopy and one 
after, in order to record the main absorption line (as in Fig. 0I c )X from which we extract the 
on-resonance optical thickness bo. The second measurement allows us to take into account the 
losses induced by the pump radiation pressure. The corresponding uncertainty on bo induces, 
at the end, a ±10% uncertainty on the critical optical thickness. Then, we fit by 



where Ai, A2, <5r and 7 are the adjustable parameters of the Raman structure described by a R , 
and the adjustable line parameterized by m,p is used to fit the background of the transmis- 
sion spectrum. This expression is not rigorous and one could instead search for the complete 
expression of the atomic response. However, it is difficult to take into account the complete 
real system, including the Zeeman degeneracy and polarization effects. Using Eq. (O allows 
us to efficiently measure the Raman parameters. As shown in Fig. [2 a), the fit is very satisfac- 
tory. For most parameters, the width 7 of the Lorentzians is larger than their shift 5r, so that 
the two Lorentzians are not separated and the Raman structure looks like a dispersion profile. 
Similarly, the corresponding scattering cross-section looks like one bell-shaped curve. The ob- 
tained widths 7 are consistent with the pump elastic scattering rate. The ratio between the gain 
amplitude A 1 and the absorption amplitude A2 is approximately constant, as expected, since it 
depends only on the Clebsch-Gordan coefficients. 

Then, the Lorentzian shape of the Raman contribution to the atomic polarizability is analyti- 
cally transformed through Kramers-Kronig relations to get 
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The atomic polarizability a R is thus fully determined. 
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Fig. 2. (a) Typical experimental spectrum (red dots) and its fit (black line) around the 
Raman resonance.The parameters obtained from the fit are A\ = 0.21 (gain amplitude), 
A 2 = 0.11 (absorption amplitude), y = 0.25r = 1.5 MHz and Sr = 0.09r = 540 kHz. 
(b) Gain and scattering cross sections, computed from Eqs. d6|7t with the Raman param- 
eters deduced from the fit. (c) Corresponding critical optical thickness. The minimum is 
bg ~ 220. This set of data corresponds to the pump parameters A = — 3.4r and £2 = 3.4r. 



However, the measurement is valid for the special polarization configuration that we have 
used, whereas for a random laser, the polarization is a priori random. To get a realistic estima- 
tion of the random laser threshold, we thus have to make an average over the polarization. We 
have checked experimentally that the coefficients Ai, Ai have a sin 2 (0) dependence with the 
relative angle 6 between the pump and the probe linear polarizations. As we have performed 
all the measurements in the optimum case (with the probe polarization perpendicular to the 
pump one), it is appropriate to multiply the measured values of Im(a R ) by 1 /2 and |a R | 2 by 
3/8 (average of sin 4 (0)). The cross-sections used to determine the random laser threshold are 
thus 

c7g/(Jo= ||a R | 2 -ilm(a R ), (6) 
o k /oq = ||a R | 2 , (7) 

where a R is experimentally determined as described above [Eqs. (0[5]l and Fig. [2 a)]. An ex- 
emple of computed cross-sections is shown in Fig. |2jb). 
Then, the critical optical thickness is easily computed from 

1 27ZCJ /ON 

rib 0cr = , (8) 

\/3a sc ag 

where the correcting tj factor (coming from the extrapolation length) writes 77 = 
(£>0crOsc +4£<7o) / (^OcrCsc + 2^0h) an d yields a second-order equation in £>ocr- The solution, 
plotted as a function of 8, is reported on Fig.|2c). As expected, the minimum is located near 
the maximum of the gain cross- section, i.e. for 8 ~ —8r. 

4. Results and discussion 



We repeat the above procedure for each couple of pumping parameters {A, £2}. The Rabi fre- 
quency of the atom-pump interaction has been calibrated by monitoring the light shift of the 
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Fig. 3. Critical optical thickness &ocr as a function of the pumping parameters A (atom- 
pump detuning) and il (Rabi frequency of the atom-pump coupling). The minimum is 
around b 0a ~ 210 - 230, for A ~ 2T and Q. ~ 2 - 3r. 



main absorption line 11351 as a function of the pump intensity, as it can be seen in Fig. [He). We 
studied only the A < part, as Raman gain is independent on the sign of A (34). Moreover, 
we have been restricted to |A| > 2 because too much radiation pressure destroys the MOT for 
|A| < 2. As the random laser will automatically start with the first frequency above threshold, 
we report in Fig.[3]the critical optical thickness defined as 

i ftx (A > n)=min[6oc r (5 J A,n)]. (9) 

8 

The minimum is around bo a ^ 210 — 230, obtained for A ~ 2r and Q. ~ 2 — 3r. 

Once the critical optical thickness is computed, the self-consistency of our model has to be 
checked on two points. Firstly, the diffusive approach leading to Eq. (Q]) requires in principle 
that the ratio L/£ sc = bo CI x <7 sc /c7o i s substantially larger than one to be justified. Nevertheless, 
it has been shown recently that this condition has not to be strictly respected, as the diffusive 
approach gives quite accurate results down to L/£ sc ~ 1 l2D . This is approximately the value 
obtained for the optimum parameters. Note also that the correction due to the extrapolation 
length (rj factor) is not negligible, as for L/£ sc ~ 1, TJ ~ 1.6. 

Secondly, we have so far only considered the Raman resonance, neglecting the influence of 
the main atomic transition at ©o, which is valid for very large detunings A ^> F. However, since 
the optimum threshold is obtained for small detuning, this is not justified. The corresponding 
one-photon transition has no gain around 5 = and then only adds scattering. This scattering 
can be decomposed into elastic and inelastic contributions. The elastic contribution will lower 
the random laser threshold, whereas the inelastic contribution, which shifts the frequency out of 
the Raman gain curve, will yield an increase of the lasing threshold. Let us examine the effect 
of the supplementary elastic scattering. It can be evaluated by 

acl = TT^ x (TT^ x ^- (10) 

The first term is the total scattering cross-section of a two-level atom, taking into account the 
detuning. The second factor, where s — 2Q. 2 /(F 2 +4A 2 ) is the pump saturation parameter, de- 
scribes the reduced scattering cross-section, keeping only the elastic part H351 . As a change of 
Zeeman sublevel is possible during a scattering event, an additional weighting factor, estimated 



as ^ ~ 0.5 l36ll , is necessary to select true elastic scattering. Adding <7 e i to a sc [Eq. (O] low- 
ers the critical optical thickness to bo CI ~ 120 — 130, with approximately the same optimum 
pumping parameters. This is however an optimistic evaluation, as inelastic scattering has not 
been taken into account. On the contrary, a conservative evaluation can be obtained by consid- 
ering inelastic scattering as pure losses, i.e. as a negative contribution to the gain cross-section. 
This is pessimistic because those photons may not be definitively lost, as further inelastic scat- 
tering can shift their frequency back on the gain curve. As previously, the inelastic scattering 
cross-section can be evaluated by 



Oinel : 



i+4A 2 /r 2 



x(l 



(11) 



The first term in the squared bracket is associated with Raman inelastic scattering whereas the 
second term is due to incoherent scattering of the two-level atom (37). Subtracting Oj ne i to the 
gain cross-section of Eq. © increases now the critical optical thickness to bo cr ~ 215 — 230. 
The optimum parameters are then located near Q. ~ 3 — 4F and A ~ 3 — 4F. Except for small 
A, where inelastic scattering is dramatic, the result is not very different from the one presented 
on Fig. [3] Especially near the optimum parameters (Q. = A ~ 3 — 4r), the optimistic evaluation 
leads to bo a - ~ 165 — 180, which is not very different from the pessimistic result (Z?ocr ~ 215 — 
230). Therefore, we conclude that the value bo a ~ 200 gives the correct order of magnitude. 
Such a high optical thickness is achievable, for instance by using compression techniques of 
magneto-optical traps l38l [39l . The corresponding ratio L/£ sc is on the order of 2. 

Random lasing occurs at a detuning from the pump \S\ ~ (5r, typically smaller than 1 MHz. 
This makes the detection of such a random laser very challenging, as the corresponding fluores- 
cence cannot easily be separated from the pump-induced fluorescence. Nevertheless, the narrow 
Raman structure could be revealed by a beat note experiment, as in l40l . or alternatively by the 
intensity correlations in the fluorescence, measured either by a homodyne technique RD or 
with a time correlator J42]. In this last experiment, a contribution from Raman scattering has 
been measured, consistent with the theoretical predictions of (36). It seems reasonable to ex- 
pect this signal to have different behaviors below and above threshold, but this remains to be 
checked by further theoretical studies. 

Finally, let us mention that our model contains several limitations, so that the numbers should 
be considered as first-order estimates. Our description of Raman gain is quite simplified in 
order to have an efficient data analysis procedure, leading to quasi-analytical results. Precise 
modelling of the complete atomic response is indeed not the goal of this article. On the contrary, 
Raman gain is used as a convenient example to illustrate the method, which is general and 
could be used with any gain mechanism, by numerically computing the real part of the atomic 
polarizability from the experimental transmission spectrum, via Kramers-Kronig relations. 

Our hypothesis of homogeneous atomic density and monochromatic and homogeneous 
pumping could also be discussed l43l . especially when high optical thicknesses are involved, as 
the pump attenuation may become important. These effects could be taken into account in nu- 
merical simulations of light transport in active and disordered medium, but have to be neglected 
to allow the analytical resolution of the diffusion equation leading to Eq. (fl} |fl~|[T3]. Note how- 
ever that the on-resonance optical thickness is not the relevant parameter for the pumping field, 
since the pump is detuned and is saturating. Moreover, diffused pump light, which penetrates 
into the sample much deeper than the coherent transmission, has also to be taken into account. 



5. Conclusion 

We have presented a method to experimentally determine the threshold of random lasing in a 
cloud of cold atoms. In this specific system, the threshold is related only to the complex atomic 



polarizability, which can be fully characterized by spectroscopic measurements. We applied 
this idea with Raman gain between light-shifted Zeeman sublevels of rubidium atoms. From 
our measurements, we estimate the critical optical thickness to be on the order of 200, which is 
achievable with current cold-atom experiments. 

The obtained critical optical thickness is lower than the one obtained with Mollow gain 
I2TI . This is in agreement with the intuition that more complex gain mechanisms offer more 
degrees of freedom, which is of course necessary to optimize several quantities (scattering and 
gain) at the same time. We are then confident that even lower thresholds can be obtained with 
other, more complex gain mechanisms, for example non-linear parametric gain induced by non- 
degenerate four-wave mixing. This may be the subject of our future investigations. 
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